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1 Introduction 

In the research on biomagnetic measurement the 
magnetic field perpendicular to the body surface 
(normal component) is widely used. However, there 
are problems of separating multiple sources 
overlapping in time when many distinct areas of 
cortex and/or cardiac tissues are active [l]-[2]. In 
these cases, the magnetic field distribution of 
normal component is not helpful for estimating the 
location and number of sources, owing to the lack 
of a dipole pattern. From the viewpoint of wave 
form analysis, it is also desired a certain method to 
detect the transient weak signals hid by the intense 
signals, such as late potential component generated 
by the fragmented activity of myocardial tissues, 
existing at the offset of the QRS complex [3]. 
Concerning to overcome the previous problems, 
new methodology is necessary to develop in both 
of not only signal processing but also SQUID 
instrumentation. There are several reports related to 
a three-dimensional (3-D) vector measurement and 
system [4]-[9]. We have developed two types of 
3-D SQUID system. One is a whole-head SQUID 
system for MEG with 195-channel and first-order 
gradiometer [10]. Other is a prototype of 3-D 
SQUID system, for MEG and MCG, with 
39-channel and second-order gradiometer [11]. 

In this paper we focus the usefulness of 3-D vector 
measurement for estimating multiple sources in 
MEG data, which are obtained by a prototype of 
SQUID system. Lastly, we introduce recent studies 
of 3-D MCG on our articles. 

2 Characteristic of a 3-D measurement for 
estimating multiple sources 

A 3-D vector measurement of MEG/MCG is 
carried out using a 3-D second-order gradiometer 
shown in Fig.l (a). Arrangement of pick up coil is 
shown in Fig.l (b) and (c). The 3-D second-order 
gradiometer is wound with Nb-Ti wire on a 
rectangular solid of 3x3x6 cm. 

The magnetic field output of this gradiometer set 
in spherical coordinate system is expressed as 
follows: 

Br = B(r,0,(p)-2B(r + Ar,0,cp)+ B(r + 2Ar,0,(p) (1) 
BO = B(r,e,<p)-2B{r,O + A0,<p) + B{r 9 0 + 2Ae,<p) (2) 


Bcp = B(r,0,cp)-2B(r,0,cp + Acp) + B(r,0,cp + 2Acp) (3) 

where position (r,0, <f) is the lowest place of the Br 
component coil shown in Fig. 1(a). Ar is a base 
line (2.8cm) of a Br coil. A0 and A(|) are a base 
line (1.4 cm) of a coil for B0 and B(|). 



(a) 3-D second-order gradiometer 


Figure 1: Configuration of a prototype of 3-D 
second-order gradiometer and arrangement of pick 
up coil for MEG and MCG. 

To discuss the magnetic field distribution (MFD) 
on the scalp measured with 3-D second-order 
gradiometer, simulation study of MEGs generating 
by multiple sources in the brain was done. MFDs, 
as shown Fig. 2 to Fig. 4, were calculated with a 
homogeneous spherical conductor model (85 mm 
radius) as a head. Fig. 2 (a), (b), and (c) show 
example of isofield contour maps of three 
components (Br, B0, B(])) for two current sources 
(each dipole moment : 30nAm) with the opposite 
direction, at D5 and F5, located at a depth of 2 cm 
from the surface of the head. Fig.2 (d), (e), and (f) 
show the 3-D MFD maps (absolute value) of three 
components. The separation of the dipoles is 4 cm. 



(b) MEG 



(c) MCG 
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This source configuration is important for MEG 
measurement since electric potentials are 
attenuated at the scalp while magnetic fields are not. 
As shown in Fig. 2(a), the map of Br significantly 
misses the dipole pattern and shows one extreme of 
-570 fT at E5. However, the map of B0 shown in 
Fig. 2(b) clearly shows two extremes of 89 fT at 
D5 and -89 fT at F5 with opposite polarity. 
Positions of D5 and F5 correspond to the place 
located each current dipole represented by arrows 
as shown in Fig. 2(a). Fig. 3 shows example of 
isofield contour maps for two current sources with 
the orthogonal direction, as shown in arrows at D4 
and F6, separating 4.2 cm apart. In this case, MFD 
of Br component (Fig.3 (a) and (d)) shows a dipole 
pattern having two extremes of 310 fT at D2 and 
-310 fT at H6 with opposite polarity due to sum of 
two current sources. This is an example of 
miscounting case as a single source to two sources 
when we calculate the inverse problem by referring 
only MFD obtained from Br component, which is 
widely used. However, on maps obtained from 
tangential magnetic fields, B0 (Fig.3 (b) and (e)) 
shows one extreme of 65 fT at F6 and B(f) (Fig.3(c) 
and (f)) also shows one extreme of 65 fT at D4. 
Positions of D4 and F6 correspond to the place 
located each current source. Namely, the number 
and/or the location of the dipoles (sources) are 
given from tangential magnetic field extremes since 
the tangential magnetic field shows maximum or 
minimum fields just above the location of a single 
dipole [8]. 

Fig. 4 shows example of isofield contour maps for 
three current sources [2] located at E3, F5, and E7. 
Each source is located 3 cm apart from E5. Br 
(Fig.4(a) and (d)) shows only one extreme of about 
750fT at E5. However, B0 (Fig.4 (b) and (e)) 
shows one extreme of about 115 fT at G5 and B(|) 
(Fig.4(c) and (f)) shows two extremes of -llOfT at 
E3 and 110 fT at E7 with opposite polarity. We can 
also predict the number and the location of sources 
by referring MFDs from B0 and B(|). 

A 3-D MEG measurement we proposed here is very 
useful for calculating inverse problem with 
multiple sources in the brain since it gives 
constraint conditions by referring not only MFDs 
of normal component (Br) but also that of 
tangential component (B0, B(|>). 

3 3-D MEG measurement and estimation 
of multiple sources overlapping in time 

First, in order to make multiple sources in the brain, 
we carried out 3-D MEG measurement of mixed 
AEF and SEF overlapping in time. The AEF was 
elicited by 1 kHz tone bursts of 50 ms duration to 
the right ear. The SEF was elicited by electric 
pulses of 0.2 ms duration with 6 to 8 mA to the 




(d) Br 


(e) B6 




Figure 2: Isofield contour maps and 3-D maps by 
two sources oriented in the opposite direction. 




Figure 3: Isofield contour maps and 3-D maps by 
two sources oriented in the orthogonal direction. 




(d) Br (e) B0 (f) Bc|) 

Figure 4: Isofield contour maps and 3-D maps by 
three sources. 




































































median nerve of the right wrist. The interval of 
stimulation was 500 ms. Stimulation to the median 
nerve was delivered 40 ms later than to the ear. 
There were two evoked fields, namely mixed AEF 
and SEF, overlapping in time. There were thirty 
measurement positions on the scalp. All magnetic 
data were averaged for 400 measurements at each 
position. The sampling interval was 0.5 ms. 
Band-pass filter was used in the range of 0.5 to 
40 Hz. 

Second, we demonstrated a multiple source 
estimation using the moving dipole inverse solution. 
Multiple sources are estimated by minimizing the 
cost function, which are the least-squares fit 
between the measured magnetic field and 
calculated magnetic field [9]. 

Fig.5(a),(b), and (c) show an example of isofield 
contour maps at 120 ms latency of AEF and at 80 
ms latency of SEF. Solid lines (positive field) of 
the Br component show the outflow of magnetic 
flux from the scalp. These lines (positive field) of 
B0 and B^) show magnetic flux directed along 


latitudes and longitudes. As shown in Fig.5(a), the 
contour map of Br shows one extreme at 
measurement position D4. It takes a note that this 
map is similar to Fig. 2(a). This map was not 
helpful in estimating the location and the number 
of sources. By referring to contour maps of B0 and 
B(|) shown in Fig.5(b) and (c), we can clearly see 
two extremes with opposite polarity on both 
contour maps. Namely, these magnetic fields 
consist of two sources with opposite polarity 
underlying each extreme [8]. We can obtain initial 
parameters of the number, the locations and the 
directions of sources by visual inspection from 
these contour maps. We carried out source 
estimation using these initial parameters. Fig.5(d), 
(e), and (f) show isofield contour maps calculated. 
The calculated maps B0 and B^) were estimated 
from the result of the inverse solution from Br. The 
goodness-of-fit of Br is 92.4%, and the 
goodness-of-fit of B0 and B^> are 88.7% and 88.6% 
respectively. 

Fig.6 shows the MRI of a subject’s head. 





(a) Br (31 OfT/In) (b) B0 (93fT/In) 


(c) B(f> (74fT/In) 



(d) Br (25OfT/In) (e) B0 (78fT/In) (f) B(f> (62fT/In) 

Figure 5: Isofield contour maps consisting of mixed AEF (120 ms) and SEF (80 ms) overlapping in time, 
(a), (b), and (c) are MEG measured, (d), (e), and (f) are MEG calculated. 






(a) Coronal 


(b) Sagittal 


(c) Transverse 


Figure 6: MRI of a subject’s head showing the superimposed source localization of AEF (open square) 
and SEF (open circle) 


Estimated sources are superimposed on the MRI. 
The localization of the AEF source (open square) at 
120 ms latency on the MRI was estimated at the 
gyrus of the sylvian sulcus with a depth of 1.8 cm 
below the scalp, which is in the auditory area. 
Localization of the SEF source (open circle) at 80 
ms latency on the MRI was estimated at the gyrus 
of the central sulcus with a depth of 1.4 cm from 
the scalp, which is the somatosensory area. The 
results of this study, analyzing magnetic fields over 
the human primary cortices, agree in general terms 
with the related findings in human studies. 

It has been applied 3-D MCG measurement to 
detection of high-frequency components of the 
QRS complex [12] and frequency analysis of 
arrhythmia [13]. 
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